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Abstract 

Engineering the optical properties using artificial nanostructured media known as metamaterials 
has led to breakthrough devices with capabilities from super-resolution imaging to invisibility. In 
this article, we review metamaterials for quantum nanophotonic applications, a recent development 
in the field. This seeks to address many challenges in the field of quantum optics using recent 
advances in nanophotonics and nanofabrication. We focus on the class of nanostructured media 
with hyperbolic dispersion that have emerged as one of the most promising metamaterials with a 
multitude of practical applications from subwavelength imaging, nanoscale waveguiding, biosensing 
to nonlinear switching. We present the various design and characterization principles of hyperbolic 
metamaterials and explain the most important property of such media: a broadband enhancement 
in the electromagnetic density of states. We review several recent experiments that have explored 
this phenomenon using spontaneous emission from dye molecules and quantum dots. We finally 
point to future applications of hyperbolic metamaterials of using the broadband enhancement in 
the spontaneous emission to construct single photon sources. 

1 Introduction 

The search for new materials naturally led to the concept of a metamaterial where the bulk electro- 
magnetic properties arise not only due to the composition but also because of the underlying structural 
resonances and near field coupling between the designed subwavelength building blocks [I]. The simul- 
taneous progress in the field of nanofabrication has been a major impetus for this unprecedented control 
over the microscopic structure. The availability of metamaterials led to exotic applications such as neg- 
ative refraction [2] , subwavelength resolution imaging in the near and far field HI [5j |6j [7j , invisibility 
devices [8j and perfect absorbers [9]. 

The class of artificial media which has emerged as one of the most important at optical frequencies 
are hyperbolic metamaterials (HMMs). These non-magnetic media have a dielectric tensor which is 
extremely anisotropic and cannot be found in nature at optical frequencies. The name is derived from 
the isofrequency curve of the medium which is hyperbolic as opposed to circular as in conventional 
media. The work in optical hyperbolic metamaterials was initiated by their ability to exhibit negative 
refraction [10] and index [H] but since then, it has diversified into a multitude of applications from 
subwavelength imaging [4j, nanoscale waveguiding and light confinement [T2J, [13], biosensing [14] to 
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enhanced nonlinearities [15]. The primary reason for this is the highest figure of merit amongst all 
other metamaterials, ease of nanofabrication, bulk three dimensional non-resonant response, deeply 
subwavelength unit cells and tunability across a broad range of wavelengths [16] . 

Historically, hyperbolic media were first explored in the field of plasma physics [17]. By applying 
a strong magnetic field to a plasma, the dielectric tensor could be made extremely anisotropic leading 
to hyperbolic isofrequency curves. In fact, negative refraction and planar focusing due to hyperbolic 
dispersion were predicted in anisotropic plasmas [T7] before Veselago's seminal paper on the negative 
index lens [18]. Other characteristic features such as the nature of dipole radiation and propagation in 
hyperbolic media were also addressed by the pioneering work of L. Felsen [19J. Inspired by the work 
in plasmas, low frequency metamaterials were designed to emulate many of the hyperbolic properties 
using microwave circuits that can be ascribed an effective index [20j . 

This article will focus on hyperbolic metamaterials at mid-infrared (mid-IR) and optical wavelength 
ranges. The range of applications at high frequencies such as subwavelength imaging or nanowaveguiding 
are fundamentally different and so too are the challenges. The purpose of this article is two-fold. First, we 
provide a comprehensive design perspective on hyperbolic metamaterials from the visible to mid-infrared 
wavelengths. We compare and explain the choice of materials and design parameters which have been 
used for various experiments. This should help experimentalists design hyperbolic metamaterials for 
specific device applications. We focus extensively on the validity of the effective medium approximation 
to describe hyperbolic metamaterials. The underlying theme of the review is the emerging area for 
plasmonics and metamaterials research: quantum applications [21]. 

It was recently predicted by Z. Jacob et al that hyperbolic metamaterials have a large photonic 
density of states leading to a broadband Purcell effect useful for applications such as efficient single 
photon sources [22l [23J [24J. This was later experimentally [25j [26], [27] and computationally [28], [29] 
verified by multiple groups, proving that hyperbolic metamaterials are a viable route for radiative decay 
engineering. The metamaterial approach of modifying spontaneous emission is fundamentally different 
from methods that involve resonant cavities [30], EI], photonic crystal structures [32], [33], slow light 
waveguides [34J or low mode volume plasmons [35j. The broadband enhancement in the spontaneous 
emission occurs due to coupling with electromagnetic states unique to the hyperbolic metamaterial. 
We review in detail the various experiments that probed the predicted effect and compare it with the 
theoretical results [TBI EE] . 

The main application of metamaterial based quantum nanophotonics would involve the coupling to 
quantum emitters such as nitrogen- vacancy centers in diamond [2H EI] . Their broadband emission are 
not ideally suited to the optical cavity based approach for enhancement. In spite of the absorption losses, 
an optimized design can lead to efficient single photon sources based on hyperbolic metamaterials [22lf38] . 
The issue of spontaneous emission near a hyperbolic metamaterial is also pertinent to the problem of loss 
compensation. A large number of emitters would emit preferentially into metamaterial modes leading to 
active photonic devices. Finally, we expect this review to pave the way for future work on non-classical 
light interaction with metamaterials. 

2 Design and Characterization of the HMM 

2.1 Isofrequency relation 

The fundamental property of the HMM is that it has dielectric properties (e > 0) in one direction but 
metallic properties (e < 0) in an orthogonal direction. The dielectric response of the metamaterial is 
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Figure 1: &-space topology. The isofrequency contour for (a) an isotropic dielectric is a sphere, and for 
(b) extraordinary waves in a uniaxial medium with extreme anisotropy (e x = e y > and e z < 0) it is 
a hyperboloid (type I HMM). (c) Hyperboloid of a type II metamaterial when two components of the 
dielectric tensor are negative (e x = e y < and e z > 0). This medium behaves like a metal reflecting all 
waves below a minimum parallel wavevector k min . The metamaterials in (b) and (c) can support waves 
with unbounded wavevectors (dashed red arrow) as opposed to an isotropic medium. 
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but with extreme anisotropy e xx • e zz < 0. Such a property does not exist in nature at optical frequencies 
and has to be engineered using nanostructured metamaterials. The word hyperbolic metamaterial arises 
from the unique isofrequency surface of extraordinary waves in such a medium which is given by 
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Note that the isofrequency surface is simply a sphere, as seen in[T] (a), for an isotropic medium (e xx = 
£yy = tzz), but opens into a hyperboloid when the dielectric constants are of opposite signs along the 
principle axes. Two possibilities exist depending on which direction is metallic leading to either a single- 
sheeted hyperboloid or a double-sheeted hyperboloid (see[j] (b) & (c)). We classify these two systems 
based on the number of negative components of the dielectric tensor as type I (e xx = e yy > while 
e zz < 0) and type II (e xx = e yy < while e zz > 0) hyperbolic metamaterials. 

The most striking feature of such metamaterials, evident from the isofrequency surface, is the exis- 
tence of unique metamaterial states which can have wavevectors far exceeding the free space wavevector 
(k = uj/c). They are referred to as the a high-fc" states of a hyperbolic metamaterial and lie at the heart 
of many device applications from imaging to quantum nanophotonics. These high- A; states are evanes- 
cent in vacuum and decay away exponentially. In stark contrast, these waves are allowed to propagate 
within an HMM. The properties and applications of these metamaterial states are discussed in Section 3 
and 4 respectively and we now focus on practical designs which can achieve the desired electromagnetic 
response. 
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Figure 2: Hyperbolic media: A materials perspective (a) Hyperbolic metamaterials can be made by 
a variety of plasmonic materials tailored to different regions of the electromagnetic spectrum from the 
visible to mid-infrared frequencies (b) multilayer realization consists of alternating subwavelength layers 
of metal and dielectric (c) HMM based on metal nanowires in a dielectric host. The spontaneous emission 
from dye molecules near these metamaterials were studied in references [271 EH] . 



2.2 Practical designs and materials perspective 

There are two main approaches to achieving the desired hyperbolic isofrequency surface using meta- 
materials. The first one relies on a metal-dielectric superlattice with subwavelength layer thicknesses 
[U [5j [39l [IB]; shown in [2] (b). The second approach utilizes metallic nanowires in a dielectric host 
[361 HU1 UU H2J H31 HI] as seen in [2] (c). Effective medium theory (see appendix A and B) shows that 
these two structures can achieve the required extreme anisotropy. It is well known that the dielectric 
constants of nanostructured metal, inevitable in hyperbolic metamaterial designs, are different from the 
established values of their bulk counterparts [45l [46]. Higher absorption in the metallic thin films and 
nanowires poses a major challenge for hyperbolic media and metamaterials in general. Deposition of 
smooth thin films of silver (for the multilayer design) with low loss can be achieved using a germanium 
wetting layer [471 148] : while the quality of thin metallic rods (for the nanowire design) has been shown 
to be significantly improved by annealing [43]. In fact, hyperbolic metamaterials have the highest figure 
of merit among any demonstrated metamaterials. 

The above mentioned designs can be tuned to function in all wavelength ranges of interest from the 
UV, visible, near-IR to mid-IR with the appropriate choice of metal and metallic filling fraction. We note 
that the key design aspects are the choice of the metal and its plasmonic response. In the visible range, 
silver forms the best choice because of its low optical losses, however as the wavelength increases we need 
to resort to other metals in the near-infrared (telecommunication wavelength) and mid-infrared regions. 
This is primarily because of the large reflectivity of the metal at higher wavelengths (lower frequencies) 
which leads to a large impedance mismatch with other media and thus reduces transmission. The best 
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Figure 3: EMT predicts that a metal dielectric composite can behave as an effective dielectric, effective 
metal, type I HMM or type II HMM depending on the wavelength and fill fraction of the metal. Optical 
phase diagrams for (a) Ag/Al 2 Os multilayer system (b) Ag/Ti0 2 multilayer system (c) AlInAs/InGaAs 
multilayer system in the mid-IR region (d) Silver nanowires in an alumina matrix. 



alternative in the long-wavelength region is a highly doped semiconductor which can behave as the 
plasmonic building block for hyperbolic metamaterials [IB] . The doping controls the plasma frequency 
of the metal and for III-V semiconductors can be tuned up to the mid-IR. At the near-IR wavelength 
ranges, transparent conducting oxides and transition-metal nitrides have emerged as an alternative 
plasmonic material to silver [321 ED]. In figure [2] (a), we outline different materials that can be used as 
building blocks for hyperbolic metamaterials. 

For the choice of dielectric, a high refractive index is required to achieve a type I response as well 
as to get a high throughput. One of the options is titanium dioxide (Ti0 2 ) in the visible spectrum 
which has a high dielectric constant. Close to UV wavelengths, however, Ti0 2 has high absorption and 
so alumina (Al 2 Os) becomes a better alternative. In figure [3j we show the optical phase diagram for 
the hyperbolic metamaterial where the effective medium response is plotted with the wavelength and 
the fill-fraction of metal. We consider both the multilayer and nanowire designs as well as visible and 
mid-IR wavelength ranges. The design and operating point for specific applications can be chosen from 
this graph. We note that all structures show both type I and type II hyperbolic metamaterial behavior 
in a broadband wavelength region. This is primarily because of the non-resonant nature of the desired 
electromagnetic response. The homogenization also predicts this metal-dielectric mixture to behave like 
an effective metal or effective dielectric where interesting properties are not expected. The transitions 
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Figure 4: Comparison of EMT and the transfer matrix (TM) method. The transmission is shown (a) 
at 440 nm, for a type I HMM consisting of a 6-layer Ag/Ti02 system (15 nm each) (b) at 750 nm, for a 
type II HMM with 0.7 fill-fraction (35nm/15nm). The insets show the reflectivity for the corresponding 
unit cells. 



between the various surfaces in &-space occur within the borders of the different coloured regions in 
the graphs [27] . The transition point where the four effective materials meet represents the plasmon 
resonance of the metal-dielectric metamaterial composite. 

2.3 How effective is effective medium theory (EMT)? 

A natural question arises whether the optical phases and hyperbolic isofrequency surfaces predicted by 
effective medium theory are achievable in the practical structure. Homogenization is only valid when the 
wavelength of light is much larger than the unit cell of the metamaterial. In hyperbolic metamaterials, 
the interesting properties arise due to propagating waves that have an effective wavelength inside the 
medium which is much smaller than the free space wavelength (high- k states). We therefore compare the 
results of effective medium theory and a practical multilayer realization of the hyperbolic metamaterial 
to validate the predictions. Throughout the review, we consider the effect of absorption, dispersion, unit 
cell size, total sample size and also the role of the substrate. Our conclusions apply for applications in 
imaging as well as quantum optics and, therefore, should guide experimentalists to choose the optimum 
design for an experiment. 

2.3.1 Propagating and evanescent wave spectrum: 

Routine characterization of metamaterials proceeds by measurement of both the phase and amplitude 
of the reflected light. Ellipsometric characterization along with the knowledge of the uniaxial nature of 
the multilayer or nanowire structure can uniquely identify the components of the dielectric tensor. We 
note that type I HMMs always have lower reflections than type II because of the positive tangential 
component of the dielectric constant. Comparing the reflected spectrum of propagating waves in the 
insets of figure |HJ we see a good match between the effective medium predictions and the practical 
multilayer structure for both type I and type II metamaterials. Irrespective of the fill-fraction, excellent 
agreement is obtained as long as the unit cell is significantly subwavelength. Surprisingly, even the unit 
cell reflectivity starts agreeing with the EMT prediction (figure 4 inset), but this is not always the case 
(figure 6 inset). 
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Figure 5: Comparison of the transmission of (p)-polarized propagating and evanescent waves through 
the metamaterial calculated using (a) transfer matrix for an 8-layer system consisting of Ag and Ti02 
(15 nm each) and (b) effective medium theory. The colour bar corresponds to the transmittance in 
log scale. The plasmon resonance occurs at a wavelength of X p = 490 nm. The multilayer structure 
achieves a type I response for A < X p and a type II response for A > X p . The existence of multiple high- A; 
metamaterial states away from the plasmon resonance is seen in the multilayer realization, in agreement 
with effective medium theory. 



Propagating waves can only probe the "low- A;" spectrum of waves allowed within the HMM. Since 
the most interesting properties stem from high-A; propagating waves within the metamaterial, it becomes 
necessary to experimentally and theoretically understand the tunneling of evanescent waves through the 
multilayer structure. In figure [4j we note the presence of multiple peaks in the transmission of large 
wavevector waves which would normally be evanescent and simply decay away in any dielectric medium. 
These correspond to coupled surface plasmon polariton waves propagating on the thin metallic surfaces 
described in the effective medium theory as high-A; states. Deviations from effective medium theory 
are natural and expected for incident evanescent waves because the effective wavelength of such waves 
starts approaching the unit cell size. Yet the key feature of high- A; modes predicted by EMT is clear in 
the practical multilayer structure, figure [5] demonstrates a similar comparison between the multilayer 
metamaterial (a) and the effective medium prediction (b) across the optical wavelength range. We note 
that the high- A; modes in both the type I and type II metamaterials converge at the plasmon resonance 
wavelength. In addition, we see a noticeable cut-off region for the type II HMM in which no light is 
transmitted through the metamaterial. Transmission is allowed between the k min point and the unit 
cell cut-off point. 

2.3.2 Number of layers: 

The effective medium limit is normally attained only in the limit of a large collection of unit cells, (N 
oc). We show that both the reflected and transmitted waves in the multilayer realization start matching 
the effective medium predictions with as little as a total of 8 layers. Starting with the unit cell, the 
addition of metal layers is seen to add more high- A; modes as predicted by EMT (see figure|6]). This is due 
to the coupled plasmonic states which increase in number on adding more metal layers. This emphasizes 
the fact that a small number of unit cells deviates significantly from the EMT prediction of the existence 
of high-A; waves. In the ideal limit, EMT predicts an infinite number of high-A; modes, because waves 
with arbitrarily large wave vectors can propagate in the hyperbolic metamaterial. However, in reality 
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Figure 6: Transmission dependence on the number of layers of a multilayer HMM (Ag/Ti0 2 15 nm 
each), (a) Type I metamaterial (A = 440 nm) for varying number of layers. Inset: Low reflectivity 
characteristic of type I HMM which has (e\\ = e x = e y > 0). (b) Type II HMM (A = 750 nm) which 
shows the presence of surface plasmon polaritons as well as high- A; modes. Inset: Reflectivity closely 
resembles that of a metal. The k m i n value of 4.2/cq agrees with EMT. Notice that the unit cell in both 
cases shows no high- A; states unlike the multilayer structures. 



the finite size of the unit cell as well as the losses curtail the effect. Specifically, we see in figure [6] that 
the optimal transmission is obtained for the 6 layer system in the type I HMM, while the type II HMM 
shows optimal performance with a minimum of 8 layers. Type II HMM has pronounced high- A; peaks 
whereas the type I HMM has a flat transmission of high-A; states which makes it suitable for imaging 
applications. 

2.3.3 Unit cell size: 

A major limitation on the highest possible wavevector which propagates in the practical realization 
is the size, a, of the unit cell. We note that for waves inside the medium which have an effective 
wavevector k e ff ~ 1/a, the metamaterial limit is no longer valid. This is in agreement with the multilayer 
simulations which show that high wavevector transmission drops off considerably for k x ~ 1/a. For such 
high-k waves, we enter the photonic crystal limit and the waves lie at the edge of the Brillouin zone of 
the ID periodic lattice. In figure [7j we see excellent agreement between two different fill-fractions that 
have the same unit cell size and therefore have approximately the same cut-off point. Also note that as 
the unit cell size increases the cut-off point decreases as well. 

2.3.4 Effect of loss: 

Even with 20 to 30 unit cells consisting of silver and titanium dioxide, the large wavevector states are 
transmitted through the metamaterial [16]. Thus the biggest limitation to device performance arises 
from the unit cell size and the losses mainly effect the propagation length of the metamaterial states. 
This is the main reason hyperbolic metamaterials have a very high figure of merit and optimum device 
performance occurs away from a resonance where the bandwidth can be large and sensitivity to losses 
is lower. 
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Figure 7: A comparison of the p-polarized transmission for different fill-fractions. Using the transfer 
matrix method, the transmittance is calculated for a Ag/Ti02 (a) 6-layer, type I HMM operating at 
440 nm, and (b) an 8-layer type II HMM with an operating wavelength of 750 nm. The fill-fractions 
seen in the legend correspond to metal/dielectric thicknesses of 16nm/24nm, 24nm/16nm, 20nm/30nm, 
30nm/20nm. The insets show the reflectance for p-polarized light. 

2.4 Regime of operation 

The wavelength and design for operation relies heavily on the application of choice. In imaging, where 
a throughput of both propagating and evanescent waves are required for high resolution, the type I 
metamaterial which has lower reflection performs best. In fact, impedance matching becomes of critical 
concern and often the regime of operation can be chosen by using e xx (X op ) = 1 in the type I region. The 
property of type II hyperbolic metamaterials to reflect propagating waves and transmit only evanescent 
waves has also been used in imaging to reduce the image distortion of subwavelength features that 
can be caused by propagating waves [39]. Near the plasmon resonance of the multilayer structure 
£ xx —> and e zz — oo, however, this is not the optimum functional wavelength due to losses and the 
large impedance mismatch with surrounding media [5T1 152] . For the nanowire metamaterial, hyperbolic 
isofrequency surfaces can be achieved with very low metal fill fraction making it exceptionally useful for 
practical applications. A considerable body of work has focused on the e zz w region of the nanowire 
array. This is also known as the L resonance since it occurs due to the coupling of the individual 
longitudinal plasmons of the metal nanowire. Applications such as ultrafast nonlinear switching as 
well as effects like additional waves all occur near this L resonance [HI [43J. For imaging, the type I 
hyperbolic behavior in the nanowire array has been demonstrated to yield excellent results both in the 
visible and near-IR (telecommunication wavelength) [53]. In section 4, we will elucidate the optimum 
regime of operation for quantum optical applications where a decreased lifetime of emitters and a high 
throughput are the desired figures of merit. 

2.5 Extraction of effective medium parameters 

It is necessary to experimentally extract the dielectric tensor of a fabricated metamaterial to compare 
with the predictions of effective medium theory. This can be performed using just the amplitude of 
the (s) and (p) polarized light and a uniaxial model of the metamaterial. The s-polarized light has a 
circular isofrequency curve and the comparison of the measured reflectivity with the predicted value 
over a range of angles can lead to determination of the tangential component (e xx = e yy ) uniquely at 
the specific wavelength of interest. The next step is to use this value and the reflected amplitude of 
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Figure 8: W-LDOS in the near field of an HMM calculated by EMT and transfer matrix methods. 
The plots above correspond to (a) an 8-layer Ag/A^O^, 15nm/15nm system calculated via the transfer 
matrix method and (b) the equivalent result predicted by EMT. The enhancement in the density of 
states in the multilayer system is curtailed by the finite unit cell size. 

p-polarized light for a range of angles to uniquely determine the third component of the dielectric tensor 
(e zz ). A least-squares fitting procedure was used in [25j [36l [54] , which gave good results. We emphasize 
that it is necessary in experiment to characterize the individual thin metallic layer and dielectric layer 
as the dielectric constants can be different than their bulk counterparts [UJHBj. A different approach 
is to study the extinction spectrum at different angles to characterize the anisotropy [15] . For nanowire 
samples, this helps to uniquely identify the e zz « wavelength which occurs at the L resonance. 

To characterize the tunneling of evanescent waves and understand the deviations from effective 
medium theory, conventional spectroscopic ellipsometry is not sufficient. For hyperbolic metamaterials 
in particular, spatial dispersion place a critical role since the characteristic high-A; metamaterial states 
have an effective wavelength in the medium comparable to the size of the unit cell. A total internal 
reflection set up can be used to measure the tunneling of evanescent waves and to characterize the 
dispersion of these unique states [55]. The deviations from EMT are significant even for propagating 
waves at wavelengths for which e zz « [43j . 

3 Photonic Density of States 

The striking property of the hyperbolic metamaterial which is key to understanding its quantum appli- 
cations is the large number of electromagnetic states supported by it. The photonic density of states 
(PDOS) similar to its electronic counterpart is the means of quantifying this. An intuitive counting 
procedure in A;-space consists of calculating the volume between the isofrequency contours at uj(k) and 
uj(k) + Auj. For vacuum, it is equivalent to the volume of an infinitesimally thin spherical shell in 
fc-space. As opposed to this, for the HMM we get a hyperboloidal shell which has infinite volume and 
in the effective medium limit leads to a broadband singularity in the density of states [22] . The physical 
reason behind this singularity is that there is no upper cut-off to the wavenumber in an ideal hyperbolic 
medium and infinitely large number of high-A; modes are supported by it. 
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Figure 9: (a) The lifetime in the near-field for a type II metamaterial with e\\ = —1.91 + 0.2H and 
e± = 57.5 + 4.4z is shown to be lower than a thin or thick layer of silver, (b) The role of the high-fc 
states in the lifetime decrease is clear from the wavevector-resolved local density of states. 



3.1 Wavevector-resolved Local Density of States (W-LDOS) 

To quantify the above argument we need to define the density of electromagnetic states accessible in the 
near field of the metamaterial and also separate the contribution of the metamaterial modes which lead 
to the enhanced DOS. We therefore define the wavevector-resolved local density of states (W-LDOS) 
p(uj, d, k) which isolates the contribution to the LDOS of various modes (specified by the corresponding 
wavevector) at a distance d in the near field of the metamaterial. The W-LDOS can be calculated in the 
near field of an effective medium slab with hyperbolic dispersion and compared with the corresponding 
multilayer realization. The W-LDOS is defined in terms of the Green's tensor for the multilayer or 
effective medium planar slab [561 157] . The explicit expression of the dyadic Green's function is given in 
appendix C. The total local density of states p(r OJ u;), at a frequency cj , is 

p{y 0) uj ) = ^Im{Tr[^(r ,r ;cj )]}. (3) 

where c is the speed of light. 

We first consider a multilayer system consisting of Ag/Al 2 Os and compare it with the predictions 
given by effective medium theory. In figure [8] (a) and (b), the W-LDOS has been calculated for these two 
systems in the near field at a distance of 30 nm. The most interesting feature evident in the W-LDOS 
is the presence of peaks corresponding to the metamaterial modes with large wavevectors. Furthermore 
they are non-resonant in nature and are present in a broad bandwidth for both type I and type II 
metamaterials. As expected, in any practical realization, the enhancement in the LDOS is curtailed by 
the finite unit cell size. In spite of this, for deep subwavelength layers which are easily achievable by 
current nanofabrication techniques, an LDOS enhancement of over 20 is easily realizable. There is no 
contribution from propagating waves to the W-LDOS of the type II metamaterial (longer wavelengths) 
as opposed to the type I metamaterial (closer to the UV). A large increase in the LDOS is also seen 
closer to the surface plasmon resonance of the individual layers. In EMT, the parallel component of the 
dielectric tensor (e x ) passes through at this wavelength. The losses are high at this resonance and it 
is better to function at a wavelength away from it. 
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Figure 10: The coupling to the high- A; modes of a Ag/TiC>2 8-layer metamaterial for (a) type I (A=440 
nm) and (b) type II (A=750 nm) metamaterial. Insets show the effect of varying intrinsic quantum yield 
of the emitter. Note that the near field decay rate is dominated by the metamaterial high- A; states. 



4 Modification of the spontaneous emission lifetime 

Fermi's golden rule immediately shows that the enhanced local density of states has a consequence on 
the spontaneous emission lifetime of quantum emitters in the near field of the HMM. These emitters 
are best represented by point dipoles that emit both propagating and evanescent waves. In free space, 
excited fluorophores relax emitting energy only into propagating waves. But near the HMM, a new 
decay route arises due to coupling of the emitted evanescent waves to the high- A; metamaterial states. 
This significantly lowers the emitter lifetime. Furthermore, the dominant mechanism of relaxation in 
the near field is the emission into the multitude of metamaterial modes and not the propagating waves 
in vacuum. 

The very presence of these large number of modes increases the decay rate causing an enhancement 
in the spontaneous emission which is known as the Purcell effect. Note that conventionally such en- 
hancement has been studied near cavity resonances but these are fundamentally limited in bandwidth. 
In hyperbolic metamaterials, the enhancement is broadband and thus they can couple to emitters which 
have a broad spectral width. Spontaneous emission lies at the heart of all light emitters from light 
emitting diodes, lasers to single photon sources. This is the main motivation behind many of the recent 
developments in engineering the spontaneous emission using hyperbolic metamaterials. 



4.1 Coupling to metamaterial modes 

In figure |9| we show the distance dependence of lifetime for a thin metal layer, the unit cell of the 
metamaterial and the multilayer hyperbolic metamaterial. The wavelength of operation is 545 nm, the 
thicknesses are 15 nm for both the metal and dielectric layers. We achieve a type II metamaterial and 
function near the transition region. In the near field, the lifetime on the multilayer structure is lowest 
in comparison to the others in agreement with effective medium theory. It should be noted that the 
lifetime decrease at close distances for the thin metal layer is due to the well known plasmonic route. 
The local density of states for the three cases shows the presence of unique metamaterial states which 
is the primary reason for the lowest lifetime. 

The coupling to the metamaterial modes is due to the evanescent waves emitted by the point emit- 
ter and occurs only in the near field. We plot the decay rate contributions due to various routes of 
relaxation for the excited emitter in figure [lOj The emission into propagating waves is reduced as the 
distance decreases while the opposite trend is noticed for the high-A; metamaterial states. The type 
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Figure 11: The effect of losses on the wavevector-resolved local density of states, p(A, d, k x ) where d=30 
nm. The area under the curve given p(A, rf, k x ) determines the overall lifetime of the dipole emitter. 
The curves correspond to W-LDOS with different optical losses for (a) a 6-layer Ag/TiC>2 system at 440 
nm (Type I HMM) and (b) an 8-layer Ag/Ti0 2 at 750 nm (Type II HMM). Higher losses broadens the 
high- k modes shown in the two figures; however, as shown in the insets, the lifetime t/t (normalized to 
vacuum) remains approximately the same. 



II metamaterial, at longer visible wavelengths, also supports metamaterial surface plasmon polaritons 
which causes a further reduction in lifetime. The insets shown in figure [10] demonstrate the lifetime 
dependence on the intrinsic quantum yield rj of the emitter. As expected, small intrinsic quantum yield 
results in reduced coupling and a higher lifetime for any given distance of the dipole emitter. 



4.2 Non-radiative decay rates 



In figure [TTJ we study the effect of losses on the local density of states. It is seen that the metamaterial 
modes which contribute to the LDOS broaden as the losses increase, corresponding to the decreased ab- 
sorption length in the metamaterial. Yet the total LDOS (integrated over all wavevectors) is unchanged 
indicating that the role of losses is detrimental mainly to the length of propagation. This also shows that 
the enhancement in the LDOS is radiative in nature since even in the zero loss limit the metamaterial 
modes lead to the same LDOS enhancement. 

We note that a closer look at the quenching of spontaneous emission is necessary to decouple the 
radiative and non-radiative decay rates. The evanescent waves emanating from a point dipole in the 
near field of a lossy structure are completely absorbed which leads to a decrease in lifetime. This is 
conventionally termed as quenching due to "lossy surface waves" [58j E7J [59j [60] . As opposed to this, 
the evanescent waves incident on a hyperbolic metamaterial are converted to propagating waves due to 
which the lifetime decrease is in fact radiative. This is analyzed in figure [12] which shows the role of 
quenching for a dipole emitter placed at d = 3 nm from various substrates (thin metal, thick metal, 
effective medium metamaterial and multilayer system). The smooth peak at large wavevectors (figure 



12 (a) and (b)) for the thin and thick metal is indicative of a route of decay other than propagating 
waves or plasmons. In the effective medium approximation, the HMM supports modes with unbounded 
wavevectors and hence radiative modes appear inspite of the quenching. In a practical realization, 
metamaterial modes with (k x > 1/a) do not exist and the multilayer realization shows quenching as 
expected. 
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Figure 12: The effect of quenching for thick and thin silver layers, a multilayer system and an effective 
medium slab. The plots show the WLDOS for an emitter in the near-field, (a) A=440 nm, d = 3 nm 
(b) A=750 nm and d = 3 nm. The smooth peak in the WLDOS at very large wavevectors is indicative 
of quenching (lossy surface waves). The quenching is reduced at a distance of d = 30 nm (c) A=440 nm 
and (d) A= 750 nm. 



4.3 Deviations from EMT 

The effective medium approximation fails when the emitter is placed at a distance much closer than the 
unit cell size. This is because the interaction between the emitter and the substrate when d « a is 
dominated by waves with k e ff ~ (l/d) » (1 /a) for which the unit cell can no longer be considered 
subwavelength. Thus the lifetime predicted by EMT is much lower as compared to a practical realization 
for very small distances. The corresponding values for the lifetime between the EMT slab and the 



multilayer metamaterial are expressed in the squared boxes of figure 12 For distances of d=30 nm, 



there is no quenching and the effective medium theory results agree well with the multilayer realization 
(figure [121(c) and (d)). 
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5 Experiments 



A number of experiments have been performed to test the nature of spontaneous emission from dye 
molecules or quantum dots near various realizations of the hyperbolic metamaterial (HMM). This paves 
the way for understanding the quantum properties of complex media such as HMMs. We explain many 
of the properties observed in experiment and comment on possible improvements. 

5.1 Multilayer and nanowire plasmonic structures 

Both the multilayer and nanowire realization of the metamaterial have been used in different wavelength 
ranges in experiment. Note that the nanowire achieves type I hyperbolic behavior whereas the multilayer 
realization leads to type II behavior in the visible range. The nanowire structure has shown up to 5 
times decrease of the emission lifetime for a Rhodamine 6G dye placed above it. Likewise, quantum dots 
and dye molecules on a Au/Al 2 Os and AgjTiOi multilayers have shown a decrease of lifetime by around 
1.5-2. This large discrepancy in experiment can be attributed to the fact that the wavelength at which 
e zz w for the nanowire structure lies within emission spectrum of the dye molecules. This leads to 
an additional enhancement due to a plasmonic resonance effect. The absorption and emission spectrum 
of the dye molecules are also altered in the vicinity of the metamaterial from that in solution. Since 
dye molecules are often suspended in a thin epoxy (dielectric) layer, light trapping causes enhanced 
absorption. Highly reflective type II metamaterials can lead to better absorption. This does not affect 
the lifetime but can affect the net output recorded at the detector [6T] , 

5.2 Intrinsic quantum yield 

The intrinsic quantum yield takes into account the intrinsic non-radiative decay of fluorophores on 
excitation. The intrinsic quantum yield is assumed to be independent of distance as long as there are 
spacer layers separating the molecules from the metal layer. The spacer layer prohibits charge transfer 
processes which cannot be accounted for by the theory outlined above. The intrinsic quantum yield 
of dye molecules is low when compared to quantum dots thus the coupling with the metamaterial is 
reduced and the change in the lifetime is not significant (see figure [l2| insets) . 

5.3 Apparent quantum yield and roughness 

The apparent quantum yield is the ratio of the light collected in the far-field by the detector compared 
to the case when the emitter is in vacuum. The experiments were performed with the excitation source 
and detector on the same side of the metamaterial. In the ideal case, since all of the emission in the 
near field goes into high-A; states which cannot reach the far-field, we expect the apparent quantum yield 
to be significantly decreased [601 [62] . Note that this is different from quenching where the lifetime and 
intensity decrease is due to the electromagnetic energy going irreversibly into heat. 

Roughness in the metamaterial structure induces major deviations in the observed apparent quantum 
yield. The metamaterial states can be out-coupled to propagating waves in vacuum due to roughness 
and non-uniformities which are inevitable in thin multilayer stacks. We note that the apparent quan- 
tum yield for a type II metamaterial which has high reflectivity can be increased through the surface 
roughness-assisted out-coupling of the metamaterial states. One important observation is that the ap- 
parent quantum yield will be very low for type I metamaterials because of the low reflectivity and high 
transmission. 
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5.4 Distributed dipole model 

To compare the experimentally measured lifetime with the semiclassical theory, the system can be 
modeled as consisting of many dipoles (oriented randomly) interacting with the metamaterial (or metal) 
substrate. These dipoles embedded inside the thin dielectric matrix (thickness D) interact with the 
substrate and part of the reflected light is collected in the far field. The variation of emission intensity 
with time can be obtained as I if) = ^2 0<d<D 1(d) exp(— T(d)t) where 1(d) is the reflected intensity in 
the far field due to a dipole placed a distance d away from the substrate and T(d) is its lifetime [25J. 

5.5 Finite sized emitters 

One limitation of the theoretical treatment is the assumption that the emitters are significantly sub- 
wavelength and hence can be represented by point electric dipoles. While this assumption works very 
well for dye molecules such as Eu+ 3 ions [58], it fails for quantum dots which can be at least 3-5 nm 
in size. In this case, the size of the emitter itself cannot be ignored and leads to deviations from the 
theoretical prediction. We note that a finite-sized emitter leads to an upper cut-off in the wavevector 
of the emitted waves. One expects this to curtail the coupling of high- A; waves and hence reduce the 
maximum possible enhancement as well as the in-coupled power [29J. But in complete contradiction, 
the opposite would occur in practice [62J. The reason for that is in the presence of losses, the finite size 
of the emitter imposes a large wavevector cut-off which actually precludes coupling to detrimental lossy 
surface waves. Thus a finite-sized emitter even in extremely close proximity decays radiatively, emitting 
a large fraction of the power into plasmonic and metamaterial modes. For careful optimization one has 
to take into account the wavevector cut-off due to three length scales which are the size of the emitter, 
size of the unit cell and the distance of the emitter from the metamaterial. 

6 Future 

6.1 Broadband Purcell Effect 

Future work will mainly be concentrated on out-coupling the metamaterial states using appropriate 
grating structures. To conclusively separate the effect of the metamaterial states from non-radiative 
paths, an experiment in a transmission mode will be needed. The non-radiative and radiative decay 
routes can then be ascertained since the light has to pass through the metamaterial and reach the 
far field before being detected. Both time resolved and time integrated (spectral) measurements are 
necessary to identify the predicted broadband Purcell effect [22J. Purcell factors of F p = 50 over a 
broad bandwidth should be attainable by current fabrication methods. One of the first applications of 
hyperbolic metamaterial based quantum nanophotonics will be single photon sources. The broadband 
response, tunability and ease of fabrication make it ideally suitable for emitters such as nitrogen vacancy 
centers in diamond. 

6.2 Active metamaterials 

The lifetime of molecules and quantum dots are significantly decreased once embedded inside the meta- 
material. This is because of the strong overlap between the metamaterial modes and the emitters. 
The emitters within the dielectric layers comprising the hyperbolic metamaterial can also provide gain 
to compensate for the losses. Active HMMs are expected to improve the device performance in all 
applications from imaging and nanowaveguiding to nonlinear optics and sensing [5H [63] . 
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7 Conclusion 



To summarize, we presented a comprehensive review of the design and characterization procedures for 
multilayer and nanowire HMMs. They are expected to be among the most promising metamaterials 
in the optical wavelength ranges with multitude of applications from imaging and sensing to quantum 
optics. Recent progress in the quantum applications of hyperbolic metamaterials and the corresponding 
experimental results were also reviewed. 
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Appendix 



A Effective medium theory for multilayer structures 

We provide here the additional details for the theoretical model of an effective medium. As mentioned, 
the extreme anisotropic response of hyperbolic metamaterials can be realized for a multilayer system 
with alternating metal and dielectric layers, where the relative permittivities are given by e m and e d 
respectively. 

To achieve homogeneity, the thickness of the layers must be much smaller than the operating wave- 
length. In practice, we require the layer thicknesses to be at least smaller than A/ 10 for EMT to be 
valid. In this regime, the uniaxial components of the dielectric tensor, eu and e± : are given by 

£\\ = ptm + (1 - p)e d (4) 

and 

pe d + (1 - p)e m 

where p is the fill-fraction of the metal in the unit cell, i.e. p = rrr for the given metal and dielectric 
thickness t m and t d . Accordingly, the plasmon resonance given by e\\ is shifted to different wavelengths 
depending on the fill-fraction of the metal. 



B Effective medium theory for cylindrical nanowires 

Following a similar approach, we can find an equivalent effective medium representation of the nanowire 
geometry given in figure [5] (c) by using the generalized Maxwell-Garnett approach. Assuming that we 
have an isotropic distribution of cylindrical nanorods, the permittivity tensor components are given by 

e_L = ptm + (1 - p)e d (6) 

and 

= (1 ~ p)e m e d + (1 - p)e 2 d 
11 (l- p )e m + (l+p)e d ' [ ) 

The axis of the cylinder nanowires is aligned along the z, or perpendicular, direction. These expressions 
are found under the assumption that p « 1 such that the local electric field across the cross-section of 
the nanowire is homogeneous [40]. The small fill- fractions used in the nanowire system generally results 
in much lower losses in the metamaterial system since less metal is used overall. 



C Dyadic Green's function formalism 

To capture the near-field effects in nanophotonic structures, we use the dyadic Green's function formal- 
ism. The Green's tensor is related to the radiation field produced by an oscillating electric dipole source 
in free space. The result, expanded as a summation of plane waves using the Weyl identity, leads to 

• /» /» 7 7 7 7 I kj~y krjQ kgpky —\— krjg j\j £ \ 

OTT J J K Z K X \ 7 7 7 2 _ p / 
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where fci = x fe{uj lc is the angular wavevector of a medium with permittivity ei. We can decompose the 
Green's foncLn tensor into its (s)- and (p)-polanzed components . and t? p , snch that 

- — ! f dk xdk y i(k x x+k vV +k z \z\) a r m \ 
where A s and A p are the tensors 

(ky ~k x k y 
k 2 x ) , (11) 


(k x k z k x k y k 2 -fk x k z (k^ + fc^) 

k x k y k 2 z kykl ^k y k z (kl + ky) ] . (12) 

Tk x k z (k x + k y ) ^k y k z (k x + k y ) (k x + k y ) 

Next, we formalize the expression of the decay rate T by following the theory developed by Ford and 
Weber [57] . This approach yields the normalized decay rate, in terms of the wavevector-resolved local 
density of states, such that 

oo 

r f 

— = 1 — 77 + 77 Re / p(oj,d,k)dk x (13) 

1 o J 



where 77 is the intrinsic quantum yield of the emitter and T is the decay rate of the emitter in free-space. 
The wavevector-resolved local density of states, normalized to free-space, is then defined as 

d ' %) = ^kl^T^ { 2^ [(1 + ^ + rP)k ' ] +Pl{1 + rP)kl } (W) 

where p\\ and p± are the parallel and perpendicular components of the dipole moment, d is the distance 
of the dipole from the planar interface, and r s,p is the reflection coefficient for (s)- and (p)- polarized 
light respectively. 

The reflection coefficients can be found by using the Fresnel equations for simple configurations, 
however, the use of the transfer matrix method is required to obtain the necessary values in multilayered 
structures. Considering a system of TV total layers, the general expression is 

( r i ) = WT 1 ™^ ( tS Q P ) (15) 

where 

N 



i»*= n^^T 1 , as) 



m = l i I, ),D ! = (t t _ l k , ) , p, = ( J idl ] . (17) 

The multilayer system consists of layers with dielectric constant and thickness d i) while k l z is the 
perpendicular wavevector component inside the ith layer. Combining all of these results, we can calculate 
W-LDOS expression in equation (14), while using numerical integration techniques to evaluate the 
normalized decay rates in equation (13). 
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